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Abstract: We show here that a purely kinetic approach to the excitation of waves by cosmic rays in the 
vicinity of a shock front leads to predict the appearance of a non-alfvenic fastly growing mode which 
is the same that was found by Bell (2004) by treating the plasma in the MHD approximation. The 
kinetic approach we present is more powerful in that it allows us to investigate different models for the 
compensation of the cosmic ray current in the background upstream plasma. 



Introduction 

There are different perspectives from which we can 
look at the role of supernova remnants (SNRs) as 
sources of the bulk of cosmic rays: from one, they 
represent the only class of sources which are en- 
ergetically viable in terms of hosting the neces- 
sary energy, provided about 10 — 20% of the ki- 
netic energy of the expanding shell is converted 
into accelerated articles. On the other hand it has 
been known for a long time [1,2] that the standard 
model needs to be pushed to its limits and beyond 
in order to explain the maximum energies which 
are observed in cosmic rays (about 10 6 GeV for 
protons). Observationally there have been recent 
claims of detection of strong magnetic fields, as 
measured from the thickness of the X-ray bright- 
ness profiles of several SNRs. These fields would 
make the acceleration process certainly easier if 
they rearrange topologically in a way to scatter the 
particles approximately at the Bohm rate. Both the 
issues of efficient acceleration and nonlinear am- 
plification of the magnetic field lead to the need to 
develop a nonlinear theory of particle acceleration 
at SNR shocks (see [3], these proceedings). 

From the theoretical point of view, the most im- 
portant aspect consists of providing a fundamental 
explanation of the mechanisms that relate cosmic 
rays as accelerated particles to the magnetic field 
amplification which has been observed. The nat- 



ural channel seems to be streaming instability, but 
other mechanisms are also possible. 

Here we discuss a kinetic approach to waves ex- 
cited by the streaming of cosmic rays upstream of 
the shock and show that there are modes which 
may be responsible for magnetic field amplifica- 
tion by a factor of order ~ 100 — 1000 for typi- 
cal parameters of a SNR. A similar conclusion was 
also obtained by using a different approach in [4] . 

Streaming instability and compensating 
currents 

In the reference frame of the shock, cosmic rays are 
approximatively stationary and roughly isotropic. 
The upstream background plasma moves with a ve- 
locity v s towards the shock and is made of pro- 
tons and electrons. The charge of cosmic cays, as- 
sumed to be all protons (positive charges) is com- 
pensated by processes which depend on the mi- 
crophysics and need to be investigated accurately. 
For instance, if n, and n e are the proton and elec- 
tron density in the background plasma respectively, 
then one can compensate the cosmic ray charge 
Ncr£ by either requiring that m + Nqr = n e , 
thereby assuming that the electrons and protons 
have different densities (and different velocity with 
respect to the shock), or we can require that pro- 
tons and electrons in the background gas have the 
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same density and velocity and that there is a fourth 
population of particles, which are cold electrons, 
at rest in the shock frame which drift together with 
cosmic rays and cancel their positive charge. The 
first avenue was adopted in [5], while the second 
was used in [6, 7]. The latter option resembles 
more closely the underlying physical approach of 
[4] where the background gas was however treated 
in the MHD approximation, and as such described 
in terms of a single density (rij = n e ). In the ref- 
erence frame in which cosmic rays are isotropic 
(namely dfcn/dfi = 0, where fcR is the distribu- 
tion function of cosmic rays) the distribution func- 
tions of the background ions, electrons and cold 
electrons are respectively 

Mp) = -46(p - miV s )5(fi + 1) 

Tl 

fe(p) = -^5(p - m e v s )5(fj, + 1) 



fe(p) 



The cosmic ray distribution can also be written as 

(1) 



fcR.(p) = ^rr-gip), 



where g(p) is a function normalized so that 
J °° dp p 2 g(p) = 1. The dispersion relation of 
waves in this composite plasma can be written as 
[8]: 
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where a runs over the particle species in the 
plasma, ui is the wave frequency corresponding to 
the wavenumber k and Q a is the relativistic gy- 
rofrequency of the particles of type a, which in 
terms of the cyclotron frequency and the Lorentz 
factor is Q = For the background plasma 

and for the cold electrons compensating the cosmic 
ray charge one has il a sa O* . 

Carrying out the usual algebraic handling of the 
different terms appearing in Eq. 2 for the four com- 
ponents of the plasma, one obtains the following 
dispersion relation, as written in the shock frame: 
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Here we have introduced the momentum p^ = 
m p Q,*/k which is the minimum momentum of 
the particles which can resonate with a wave with 
wavenumber k. The momentum po is the minimum 
momentum of the accelerated particles. The above 
dispersion relation is considerably simplified in the 
limit of frequency u> small compared with f2* . In 
particular we shall assume that ui + kv s <C ft* <C 
|fi*|. Moreover we will neglect the displacement 
current (the unity), so that the dispersion equation 
reads 



v\k 2 = (Lo + kv s ) 2 ±un*^-{i± 
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The dispersion relation in the frame of the up- 
stream plasma is obtained by the change uj — > 
uj — kv s . Moreover, neglecting again uj with re- 
spect to kv s we get 



2,2 2 kv s n*N C R 
v A k = uj =F 



[l±h(k)±ih(k)}. 

(2) 

This dispersion relation is basically the same as 
that obtained in [4] and has therefore the same im- 
plications. The terms proportional to I 2 are those 
describing the resonant interaction of the waves 
and particles. The non resonant part is proportional 
to 1 + h. 

In the following we investigate these implications 
in the case of a power law spectrum of cosmic 
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rays g(p) oc (p/po)~ 4 - In order to have the cor- 
rect normalization of g(p) we must require g(p) = 
(l/pl)(p/po)~ 4 - It follows that the energy density 
in the form of accelerated particles is 



U C R 



CR 



. dpp 3 cg(p) = ^^-podnR 

Z Jpo Z 

(3) 

where R = p m ax/Po> an d Pmax is the maximum 
momentum of the accelerated particles. The cos- 
mic ray number density can then written in terms 
of energy density as 



N C R = 



2Ucr 
PqcAylR 



(4) 



The constant in front of the square brackets in Eq. 
2 can therefore be written as 
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where £ = 2r](v s /c)(l/ lni?) and r\ = 
Ucr/ {pu 2 ) is the fraction of the kinetic energy 
that is transformed into cosmic rays at the shock. 
rL,o is the gyration radius of protons with momen- 
tum po- 

The solution for the imaginary part of the fre- 
quency from the dispersion relation Eq. 2 can be 
easily found to be: 

wi = Im[uj] = - i (v 2 A k 2 ± ct(1 ± 7i)) ± 
y(v 2 A k 2 ±a(l±h)f+a 2 I 2 

In Eq. 2, for k <~ 1/tl,o, the two functions I\{k) 
and h{k) are of order unity. Therefore it is inter- 
esting to investigate what happens in this situation 
when a -c k 2 v\ and a » k 2 v\. 

When a k 2 v 2 A we can expand Eq. 6 in the ratio 
a/{v 2 A k 2 ) < 1 and we get: 
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This is the standard result already obtained in 
previous literature, which implies the (resonant) 
growth of Alfven waves weakly driven by the 
streaming of cosmic rays. It is however instructive 
to see what the condition a <C k 2 v A corresponds 
to for typical parameters of a supernova. Let us 
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Figure 1 : Imaginary part (solid curve) and real part 
(dashed curve) of the frequency. 
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Figure 2: Phase speed of the waves compared with 
the Alfven speed and the sound speed. 



adopt B Q = IfiG, v s = 5000 km/s, n = 1 cm~ 3 , 
Pmax — 10 6 GeV/c, po = 1 GeV/c. The condi- 
tion therefore leads to i] <C 10 ~ 6 . This is a very 
stringent condition, requiring that only a millionth 
of the SNR energy is converted into cosmic rays. 
If this were the case, SNRs would have a negli- 
gible role as sources of galactic cosmic rays. It is 
possible to change the values of the parameters and 
bring the limit on rj up to <~ 10~ 3 — 10~ 4 , but still it 
appears that the condition necessary for the growth 
to be considered as weakly driven by the cosmic 
ray streaming is too stringent. In other words, the 
request that 77 ~ 0.1 typically formulated, naturally 
leads to the conclusion that the wave growth takes 
place in the opposite regime, 3> k 2 v A . 
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In this case it is straightforward to show that u>j w 
y/B. 

There are two interesting aspects of this growing 
mode that deserve some more attention. As al- 
ready pointed out in [4], the fact that a » k 2 v\ 
implies that the mode does not resemble an Alfven 
wave growing under the effect of cosmic rays. An- 
other point is that the growth takes place also due 
to the non-resonant interaction of the wave with the 
streaming cosmic rays. 

A more detailed analysis of the modes requires to 
keep track of the wavenumber k. The imaginary 
and real parts of the frequency are shown in Fig. 
1 for B = 3fiG, V s = 10 9 cm s" 1 , 77 = 0.1, 
rij = 1cm -3 . On the x-axis we plotted fcr^o, 
while the y-axis is in units of v 2 /(crL,o)- The plot 
is limited to the most interesting region fcr^.o > 1- 
The peak in /m[w] corresponds to the maximum 
growth rate. Note that the highest k part of the plot 
corresponds to the case of resonant scattering, and 
in fact for lower values of the shock speed (already 
for v s = 10 8 cm s -1 ) this is the only mode left. 
This suggests that at different stages in the evolu- 
tion of the remnant the non-resonant or the reso- 
nant growth will be relevant. Note also that the 
maximum growth occurs at k S> l/r^o, namely 
on scales much smaller than the Larmor radius of 
the particles. The effect of these magnetic oscil- 
lations on the motion of the particles need to be 
carefully investigated. 

One final note is on the phase velocity of the 
waves: in Fig. 2 we show the phase velocity 
Re\uj]/k as compared with the Alfven speed and 
the sound speed (for a temperature of the upstream 
gas T — \Q 5 K). We can see that there is a range 
of values of k (not the ones where the maximum 
growth occurs) for which the phase speed is larger 
than the Alfven speed and for small fc's even larger 
than the sound speed. These supersonic modes are 
likely to develop shocks in the precursor that might 
inhibit the acceleration process, also due to the fact 
that the effective speed of the scattering centers, 
v s — Re[w]/k becomes small enough that the ac- 
celeration may be substantially suppressed. 



Conclusions 

We showed that the growing mode found in [4] also 
results from a kinetic approach. The dispersion re- 
lation leading to the identification of the resonant 
and non-resonant modes is basically the same in 
the two approaches. It is also shown that for typ- 
ical parameters of a supernova the magnetic field 
is likely to be substantially amplified in the shock 
vicinity as due to non-resonant growth of a mode 
that does not resemble an Alfven wave weakly 
modified by cosmic ray streaming. The growth 
of this mode is faster than the standard resonant 
growth obtained in the limit of small efficiency of 
acceleration and may make particle acceleration to 
high energies easier. 
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